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CHAPTER  I 


INTRODUCTION 

The  emphasis  in  recent  years  concerning  alternate  sources  of  energy 
has  resulted  in  a  large  body  of  research  concerning  the  solar  radiation 
received  at  the  earth's  surface.  This  surface  solar  radiation,  termed 
insolation  or  hemispheric  radiation,  is  comprised  of  two  components;  the 
direct  and  diffuse.  The  hemispheric  radiation,  variously  called  the 
total  or  total  hemispheric,  is  the  total  irradiance  received  from  the 
upward  hemisphere.  The  direct,  variously  termed  direct  beam  or  simply 
beam,  component  is  that  irradiance  received  directly  from  the  sun  (usu¬ 
ally  within  a  field  of  view  of  about  5°).  The  direct  irradiance  multi¬ 
plied  by  the  cosine  of  the  solar  zenith  angle  gives  the  contribution  to 
the  hemispheric  irradiance  from  the  direct  component.  The  remainder  of 
the  hemispheric  irradiance,  which  comes  from  all  angles  in  the  upward 
hemisphere  not  contained  in  the  direct  measurement,  is  called  the  dif¬ 
fuse  component.  One  final  definition  is  required.  The  extraterrestrial 
solar  radiation  is  that  irradiance  available  at  the  top  of  the  atmos¬ 
phere  on  a  surface  normal  to  the  local  vertical.  The  value  of  this 
extraterrestrial  irradiance  is  affected  by  variations  in  the  sun-earth 
geometry.  These  effects  include  an  annual  variation  of  about  ±3.3%  in 
the  sun-earth  distance  due  to  the  ellipticity  of  the  earth's  orbit,  the 
inclination  of  the  earth's  axis,  and  the  rotation  of  the  earth  about 
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its  axis.  These  effects  are  predictable  and  a  method  of  calculating 
the  extraterrestrial  solar  radiation  is  detailed  in  Chapter  II. 

Liu  and  Jordan  (1960)  pioneered  research  in  determining  the  rela¬ 
tionship  between  the  hemispheric  and  the  diffuse  component.  Their  work 
was  based  primarily  on  data  measured  at  Blue  Hill,  Massachusetts.  Since 
this  initial  study,  many  authors  have  investigated  these  relationships 
for  various  other  locations.  Most  notable  among  these  studies  are  those 
by  Ruth  and  Chant  (1975)  for  selected  Canadian  locations,  Orgill  and 
Hollands  (1976)  for  Toronto,  Canada,  and  Collares-Pereira  and  Rabl 
(1978)  for  locations  in  New  Mexico,  Texas,  California,  and  Massachu¬ 
setts.  Chapter  III  of  this  paper  will  deal  with  the  diffuse-hemispheric 
relationship  based  on  data  measured  at  Fort  Collins,  Colorado  --  a  high 
altitude,  dry  environment. 


CHAPTER  II 


CALCULATION  OF  EXTRATERRESTRIAL  SOLAR  RADIATION 

Normalization  of  the  hemispheric  irradiance  and  its  components  is 
accomplished  by  forming  the  dimensionless  ratio  of  the  hemispheric 
and/or  a  given  component  to  an  appropriate  value  of  the  extraterrestrial 
solar  radiation  (ETR).  It  should  be  pointed  out,  however,  that  the 
direct  component  should  be  normalized  by  that  irradiance  available  at 
the  top  of  the  atmosphere  on  a  surface  normal  to  the  sun's  rays.  For 
the  purposes  of  this  study,  this  modified  ETR  is  termed  the  direct  ETR 
(DIRETR).  The  purpose  of  this  normalization  procedure  is  to  eliminate 
the  dependence  of  the  surface  values  on  solar  position.  This  is  espe¬ 
cially  important  whenever  comparisons  are  to  be  made  for  different  times 
(i.e.  one  hour  to  another,  one  day  to  another,  etc.)  or  for  different 
locations. 

Sellers  (1965)  gives  procedures  for  calculating  instantaneous  and 
daily  extraterrestrial  solar  radiation  values  and  a  graphical  method  for 
determining  daily  and  monthly  values.  Guttman  and  Mathews  (1979)  out¬ 
line  a  method  for  calculating  values  for  variable  time  periods  of  an 
hour  up  to  one  day.  The  remainder  of  this  chapter  will  deal  with  a 
method  of  calculating  extraterrestrial  solar  radiation  over  various  time 
spans  and  with  known  errors. 
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A.  The  Method 

Instantaneous  extraterrestrial  solar  radiation  can  be  calculated 
from  the  relationship  (Sellers,  1965): 

Q  =  S(R/R)2  cos  Z  [2.1] 

where, 

Q  is  the  extraterrestrial  solar  radiation, 

5  is  the  solar  constant,  taken  to  be  1377  Watts  per 

square  meter  (U.S.  Department  of  Commerce,  1978), 

(R/R)  is  the  ratio  of  the  mean  to  actual  sun-earth  distance, 
cos  Z  is  the  cosine  of  the  solar  zenith  angle. 

Determination  of  the  relative  sun-earth  distance  parameter  ((R/R)2)  and 
the  solar  position  using  Fourier  Series  was  first  proposed  by  Spencer 
(1971).  Thompson  (1981)  devised  more  accurate  Fourier  Series  represen¬ 
tations  for  these  parameters.  The  Fourier  Series  presented  in  this 
paper  are  based  on  the  formulations  of  Thompson  (1981).  They  have  been 
modified,  however,  to  facilitate  computer  programming.  All  series 
formulations  were  based  on  data  obtained  from  the  American  Ephemeris  and 
Nautical  Almanac  for  1965.  The  Fourier  Series  formulation  for  calcu¬ 
lating  the  relative  sun-earth  distance  parameter  is  detailed  in  section 
II  A.l.  The  solar  position  is  determined  from  the  geometric  formula: 

cos  Z  =  sin  4>  sin  6  +  cos  $  cos  6  cos  h  [2.2] 

where, 

4>  is  the  latitude  of  the  observation, 

6  is  the  solar  declination,  +  for  north  latitudes, 

-  for  south  latitudes, 
h  is  the  local  hour  angle  of  the  sun. 
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The  Fourier  Series  formulation  for  calculating  solar  declination  is 
detailed  in  section  II  A. 2.  The  local  hour  angle  is  related  to  Green¬ 
wich  mean  time  (GMT)  as  follows: 

h  =  |2  (GMT  -  12  +  T)  -  «|»  [2.3] 

where, 

tj/  is  the  longitude  of  the  observation  in  radians,  +  for 
west  longitude,  -  for  east  longitude, 

T  is  the  equation  of  time. 

The  equation  of  time  corrects  for  the  difference  between  true  solar  time 
and  mean  time  which  is  in  common  use.  The  Fourier  Series  used  to  calcu¬ 
late  the  equation  of  time  is  detailed  in  section  II  A. 3. 

The  Fourier  series  formulations  are  approximations  to  the  exact 
values.  Data  was  obtained  from  the  American  Ephemeris  and  Nautical 
Almanac  for  the  years  1974,  1979,  1980  and  1981  and  compared  to  values 
calculated  using  these  series  representations.  The  errors  thus  obtained 
are  given  at  the  end  of  the  section  detailing  each  series. 

1.  Relative  sun-earth  distance 

The  Fourier  series  representation  for  the  relative  sun-earth  dis¬ 
tance  parameter  is: 

(R/R)2  =  1.000145  +  0.033382  cos  a  +  0.001952  sin  a 

+  0.000698  cos  (2a)  [2.4] 

The  time  parameter  (a)  is  given  by: 

a  =  0.0172017221  [ET  +  (N00N/24)]  [2.5] 

where, 

ET  is  the  number  of  days  from  a  base  date  of  0  January  1974 
(i.e.  1  January  1974  is  day  one), 
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NOON  is  local  noon  (in  Greenwich  mean  time)  for  the  longitude 
of  the  point  for  which  the  calculation  is  being  made, 
the  constant  (0.0172017221)  is  2n  divided  by  the  annual 
period  for  the  earth  to  return  to  the  same  distance 
from  the  sun  (365.2649  days). 

The  maximum  error  in  the  calculated  square  of  the  relative  sun- 
earth  distance  using  this  formula  is  ±  0.011%. 

2.  Solar  decl ination 

The  Fourier  series  representation  of  the  sun's  declination  is: 

6  =  0.00666  -  0.400866  cos  p  +  0.06523  sin  p 

-  0.00667  cos  (2p)  +  0.00057  sin  (2p) 

-  0.00274  cos  (3p)  +  0.00124  sin  (3p)  [2.6] 

The  time  parameter  (p)  is  given  by: 

p  =  0.0172027912  [ET  +  (TIME/24)]  [2.7] 

where, 

TIME  is  the  mid-point  of  the  integration  period  (in 
Greenwich  mean  time), 

the  constant  is  derived  from  an  annual  period  of 
365.2422  days. 

The  maximum  error  in  the  calculated  declination  using  this  formula 
is  ±  1  minute. 

3.  Equation  of  time 

The  Fourier  series  for  calculating  the  equation  of  time  is: 

T  =  0.0098  cos  p  -  0.1227  sin  p  -  0.05^0  cos  (2p) 

-  0.1581  sin  (2p)  -  0.0056  sin  (3p) 

-  0.0021  cos  (4p)  -  0.0030  sin  (4p)  [2.8] 

where  the  time  parameter  (p)  is  calculated  by  equation  [2.7]. 
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The  maximum  absolute  error  in  the  calculated  equation  of  time  using 
this  formula  is  ±  12  seconds. 

B.  The  Computational  Procedure 

Integration  of  equation  [2.1]  over  some  time  interval  for  which  the 
declination  and  relative  sun-earth  distance  can  be  assumed  constant 
yields: 

=  S(R/R)2  [sin  4>  sin  6  (h2  -  h^) 

+  cos  $  cos  6  (sin  h2  -  sin  h-^/ANGVEL  [2.9] 

where, 

ANGVEL  is  the  earth's  angular  velocity  in  radians  per  hour. 
Because  it  changes  so  slowly  the  relative  sun-earth  distance  was  cal¬ 
culated  for  noon  of  each  day  and  assumed  constant  for  that  day.  Thomp¬ 
son  (1981)  points  out  that  calculation  of  the  solar  declination  once 
every  four  hours  gives  satisfactory  results;  however,  for  this  study 
values  of  both  declination  and  the  equation  of  time  were  calculated  at 
the  midpoint  of  each  integration  interval  and  assumed  constant  for  that 
interval.  A  one  hour  time  interval  was  chosen  for  use  in  this  study, 
except  that  when  h^  was  less  than  the  sunrise  hour  angle  it  was  set 
equal  to  the  sunrise  hour  angle  and,  similarly,  when  h2  was  greater  than 
the  sunset  hour  angle  is  was  set  equal  to  the  sunset  hour  angle. 

The  appendix  is  a  program  listing  (FORTRAN  V)  of  this  computational 
scheme  with  the  time  interval  (TIMINC)  set  to  one  hour.  The  program 
will  output  daily  total  extraterrestrial  solar  radiation  for  the  years 
and  location  specified.  Values  other  than  daily  totals  (i.e.  instantan¬ 
eous,  hourly,  etc.)  can  be  obtained  by  changing  the  beginning  (BEGGMT) 
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and  ending  (ENDGMT)  times.  For  mean  monthly  values  a  simple  averaging 
loop  should  be  added. 


C.  Accuracy  of  the  Method 

Given  the  errors  inherent  in  calculating  the  relative  sun-earth 
distance,  declination,  and  equation  of  time,  the  maximum  error  in  cal¬ 
culated  extraterrestrial  solar  radiation  has  been  determined.  The 
errors  were  maximized  with  the  declination  set  at  23°  and  the  hour 
angles  (which  are  a  function  of  the  equation  of  time)  set  at  90°for  the 
instantaneous  error  determination  and  at  0°  and  15°  for  the  hourly 
determination.  Hourly  errors  were  then  summed  for  15  hourly  calcula¬ 
tions  to  arrive  at  the  daily  errors.  It  should  be  noted  that  the  im¬ 
proved  performance  of  the  daily  calculation  over  the  hourly  calculation 
is  due  to  the  fact  that  the  equation  of  time  correction  is  symmetrical 
about  true  solar  noon.  For  this  reason  whenever  the  calculation  is 
performed  over  an  interval  which  is  also  symmetrical  about  true  solar 
noon,  the  equation  of  time  correction  is  cancelled  and  the  error  associ¬ 
ated  with  it  makes  no  contribution  to  the  calculated  extraterrestrial 
solar  radiation.  It  should  be  emphasized  that  the  stated  errors  are  due 
solely  to  the  computational  procedure  and  do  not  account  for  any  changes 
in  solar  output  of  long  term  changes  in  sun-earth  geometry. 

1.  Instantaneous  values 

-2 

Maximum  absolute  error  1.3  Wm 
Maximum  percent  error  0.10% 

2.  Hourly  values 

-2  -1 

Maximum  absolute  error  4  W  hr  m  hr 
Maximum  percent  error  0.32% 
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3.  Daily  values 

-2 

Maximum  absolute  error  5  W  hr  m  day 
Maximum  percent  error  0.05% 


CHAPTER  III 


ANALYSIS  OF  FORT  COLLINS  SOLAR  RADIATION  DATA 

Both  the  direct  component  (DIR)  and  hemispheric  irradiance  (HEM) 
have  been  measured  at  the  Atmospheric  Science  Department  of  Colorado 
State  University  since  the  spring  of  1975.  The  hemispheric  irradiance 
is  measured  by  an  Eppley  Precision  Pyranometer  in  the  spectral  interval 
0.3  pm  to  3.0  pm,  while  an  Eppley  Pyrhel iometer  measures  the  direct 
component  in  the  same  wave  length  band.  Cox  and  McKee  (1978)  and  Cox 
and  McKee  (1980)  give  a  complete  description  of  the  facility  and  the 
instrumentation,  as  well  as  hourly,  daily,  and  monthly  values  of  the 
direct  component  and  hemispheric  irradiance.  Because  of  the  complete¬ 
ness  of  the  data,  the  years  1977  through  1980  were  chosen  as  the  data 
base  for  this  study.  The  diffuse  component  (DIF)  was  calculated  from 
the  relation: 

DIF  =  HEM  -  (DIRHcos  Z).  [3.1] 

The  method  described  in  Chapter  II  for  determining  solar  position  was 
used  in  calculating  cos  Z. 


A.  Data  Handling 

Quality  of  solar  radiation  data  can  be  affected  by  the  calibration 
of  the  pryanometer  or  pyrhel iometer ,  mechanical  or  electronic  problems 
with  the  recording  system,  misalignment  of  the  pyrhel iometer,  and/or 
dirty  or  snow  covered  optical  surfaces.  Therefore,  the  following 
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decision  criteria  were  applied  to  the  data  prior  to  further  analysis  and 
resulted  in  a  total  data  set  of  707  days. 

1.  Missing  data 

Data  for  an  entire  day  was  disregarded  if  for  one  or  more  hours 
either  the  direct  or  hemispheric  values  were  missing.  A  total  of  506 
days  fell  into  this  category. 

2.  Suspicious  data 

If  the  calculated  diffuse  component  was  less  than  zero  the  entire 
day's  data  was  disregarded.  This  condition  results  from  the  direct 
component  being  too  large  or  the  hemispheric  being  too  small,  either  of 
which  indicate  instrument  problems.  Only  nine  days  fell  into  this 
category. 

3.  Problems  at  sunrise/sunset 

Near  the  times  of  sunrise  and  sunset  it  is  possible  for  the  solar 
zenith  angle  to  exceed  ±90°,  in  which  case  the  cosine  of  the  solar 
zenith  angle  will  be  negative  and  the  resulting  calculated  diffuse  com¬ 
ponent  will  be  greater  than  the  hemispheric  (assuming  the  direct  is 
non-zero).  For  these  events  the  hourly  diffuse  component  was  set  equal 
to  the  hemispheric  value. 

4.  Excessively  large  values  of  diffuse 

All  hours  for  which  the  inequality 

DIF  S  0. 5(S  x  cos  Z)  [3.1] 

was  not  true  were  eliminated  (Collares-Pereira  and  Rabl ,  1978).  This 
condition  results  from  the  direct  component  being  too  small  with  respect 
to  the  hemispheric,  the  most  probable  cause  being  misalignment  of  the 


pyrhel iometer.  A  total  of  239  days  were  eliminated  for  not  satisfying 
the  above  inequality. 


B.  Results  of  Analysis 

1.  Frequency  distributions 

Table  1  gives  the  number  of  days,  by  month  for  which  the  diffuse 
component  to  hemispheric  (DIF/HEM)  ratio  fell  within  the  various  inter¬ 
vals.  Table  2  and  Table  3  are  similar,  but  for  the  diffuse  to  extrater¬ 
restrial  (OIF/ETR)  and  hemispheric  to  extraterrestrial  (HEM/ETR)  ratios, 
respectively.  Inspection  of  Table  1  reveals  a  maximum  occurrence  of 
days  (23%)  in  the  0.10-0.19  range.  These  represent  days  of  relatively 
clear  conditions.  A  secondary  peak  (11%)  is  observed  in  the  range 
0.90-1.00,  representing  overcast  conditions  when  the  daily  direct  compo¬ 
nent  is  equal  to  or  very  nearly  zero.  From  Table  2  we  find  that  on  75% 
of  the  days  the  DIF/ETR  ratio  falls  in  the  0.10-0.29  range.  Table  3  in¬ 
dicates  that  on  29%  of  the  days  the  HEM/ETR  ratio  falls  in  the  0.60-0.69 
range,  and  no  evidence  of  a  secondary  maximum  exists.  Table  4  gives  the 
frequency  distribution  of  the  direct  to  direct  extraterrestrial  ratio 
(DIR/DIRETR).  Most  notable  here  is  the  relatively  uniform  distribution 
in  the  range  of  0.00  to  0.59,  with  the  difference  in  the  maximum  and 
minimum  occurrences  being  only  about  9.8  percent. 

2.  An  extended  period  of  low  surface  solar  radiation 

The  Front  Range  of  the  Colorado  Rocky  Mountains  has  the  reputation 
of  being  generally  sunny  and  mild  and,  therefore,  an  excellent  environ¬ 
ment  for  solar  energy  applications.  The  question  arises  as  to  how  well 
deserved  this  reputation  may  be.  In  an  attempt  to  at  least  partially 
answer  this  question,  the  data  was  scanned  to  determine  the  number  of 
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times  when  the  daily  hemispheric  to  extraterrestrial  ratio  was  less  than 
0.35  for  a  period  of  three  or  more  consecutive  days.  The  criteria 
(HEM/ETR)  <  0.35  was  chosen  because  Orgill  and  Hollands  (1976)  point  out 
that  values  in  this  interval  represent  extremely  cloudy  days  with  over 
90%  of  the  total  insolation  being  diffuse.  Six  such  cases  were  found. 
They  were  27-29  April  1977,  16-18  August  1977,  6-9  February  1978, 

1-3  May  1978,  6-8  May  1980.  Because  of  the  duration  of  the  February 
1978  event,  the  fact  that  this  is  the  only  period  which  was  found  in  the 
primary  heating  season,  and  having  noted  that  the  HEM/ETR  values  for  11 
and  12  February  were  also  below  the  threshold,  a  closer  examination  of 
this  period  was  made. 

a.  The  meteorological  situation.  A  high  pressure  cell  was  cen¬ 
tered  over  southwestern  Colorado  on  5  February  with  a  stationary  front 
lying  just  east  of  Fort  Collins.  Fort  Collins  experienced  light  (less 
than  5  knots)  south  easterly  flow  throughout  the  entire  day  and  by  1900 
LST  began  observing  the  classic  upslope  conditions  (ceiling  at  100  feet 
and  two  miles  visibility  in  fog).  On  6  February  the  front  moved  west 
over  the  Rocky  Mountains.  Ceilings  were  generally  at  500  feet,  visi¬ 
bility  three  to  four  miles  in  fog,  and  winds  remained  light  out  of  the 
south  through  southeast.  A  cold  front,  having  moved  off  the  Pacific 
Ocean,  merged  with  the  stationary  front  on  the  8th  and  produced,  on  the 
7th  and  8th  of  February,  ceilings  100  to  300  feet,  visibility  h  to 
k  mile  in  fog  and  occasional  light  snow,  and  winds  calm  to  light  south¬ 
easterly.  On  the  morning  of  the  9th,  a  second  frontal  system  moved  on 
shore  and  stagnated  over  central  Nevada  and  central  California  and 
remained  in  this  location  through  the  10th  as  the  low  pressure  cell 
centered  over  central  Nevada  intensified.  The  system  began  a  slow 
eastward  movement  on  11  February  which  continued  until  around  midnight 
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when  it  merged  with  the  stationary  front  over  the  Colorado  Rockies. 
During  this  entire  period  Fort  Collins  continued  to  experience  up-slope 
conditions,  with  ceilings  generally  300  to  800  feet  (except  5,000  to 
7,000  on  the  10th),  visibilities  from  less  than  one  up  to  five  miles  and 
light  southeasterly  flow.  The  active  portion  of  the  cold  front  then 
dropped  south  into  east  central  New  Mexico,  and  by  0800  LST  the  station¬ 
ary  front  had  dissipated  and  the  surface  winds  became  light  out  of  the 
northeast,  although  the  ceiling  remained  at  between  2 00  and  500  feet 
throughout  most  of  the  daylight  hours.  By  early  morning  on  the  13th  the 
clearing  trend  was  well  established. 

b.  Radiative  observations.  The  hemispheric/extraterrestrial  ratio 
and  direct/direct  extraterrestrial  for  each  day  during  this  period  is 
given  below. 


DAY 

HEM/ETR 

DIR/DIRETR 

6 

0.25 

0.00 

7 

0.21 

0.00 

8 

0.11 

0.00 

9 

0.29 

0.00 

10 

0.43 

0.09 

11 

0.15 

0.00 

12 

0.23 

0.11 

In  reviewing  the  weather  records  for  this  February  1978  case,  it 
appeared  that  there  were  other  periods  of  extreme  cloudiness  that  had 
not  been  detected  by  the  proposed  threshold.  Time  series  plots  of 
average  cloudiness  during  daylight  hours,  the  DIR/DIRETR  ratios,  and  the 
HEM/ETR  ratios  for  the  period  5  January  to  5  March  were  made,  and  are 
given  in  Figure  1.  Three  periods  stand  out  in  the  cloudiness  and  DIR/ 
DIRETR  plots  --  12-20  January,  5-16  February,  and  24  February  through 
3  March.  Average  cloudiness  during  the  5-16  February  event  was  9.6/10, 
while  the  climatological  mean  for  February  is  only  5/10.  The  average 
DIR/DIRETR  for  the  period  was  only  0.04.  Temperature  deviations 


Figure  1.  Time  series  plots  of  average  daylight  hours  cloudiness,  average  daily 

direct/direct  extraterrestrial  ratios,  and  the  average  daily  hemispheric/ 
extraterrestrial  ratios  for  5  January  to  5  March  1978.  Points  indicated 
by  an  M  are  missing,  those  indicated  by  an  E  are  estimated. 
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(observed  average  minus  climatological  mean)  were: 

Maximum  (°F)  -13.8 

Minimum  (°F)  +2.3 

The  synoptic  situation  during  the  January  and  February/March  cases 

proved  to  be  practically  identical  to  that  detailed  for  5-12  February. 

The  results  of  further  investigation  are  as  follows: 


12-20  January 

24  February  -  3  March 

Mean  Observed 

Mean  Observed 

Cloudiness 

4.5/10 

8.9/10 

5/10 

9.9/10 

DIR/DIRETR 

-- 

0.03 

— 

0.04 

Maximum  Temperature  (°F) 

40.6 

33.0 

44.8 

36.5 

Minimum  Temperature  (°F) 

14.6 

10.0 

18.0 

20.2 

We  therefore  have  three  events  covering  29  days  in  a  total  period  of 
only  58  days  where  the  average  daylight  hours  cloudiness  exceeded  88% 
cloud  cover,  the  average  daily  DIR/DIRETR  ratio  was  only  0.04,  and 
temperatures  where  highly  variable  with  maximums  from  7.6  to  13.8  de¬ 
grees  below  the  means  and  minimums  from  2.3  degrees  above  to  4.6  degrees 
below  mean  values. 

A  four  year  data  base  is  insufficient  to  establish  a  real  solar 
radiation  climatology.  In  general,  however,  it  would  seem  that  the 
reputation  of  the  Front  Range  is  well  deserved,  though  occasional  ex¬ 
tended  periods  (three  days  or  more)  of  low  insolation  may  be  expected. 
The  January-March  1978  event  is  certainly  an  anomally  and  should  not  be 
considered  representative  of  a  normal  year. 

C.  Relationships  Between  Solar  Radiation  Components 

Subsequent  to  the  initial  work  by  Liu  and  Jordan  (1960)  relating 
the  diffuse  and  hemispheric  irradiance,  many  authors  have  investigated 
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this  relationship  for  various  other  locations.  Among  these  Ruth  and 
Chant  (1975),  Orgill  and  Hollands  (1976),  and  Collares-Pereira  and  Rabl 
(1978)  found  correlations  which  deviated  significantly  from  those  of  Liu 
and  Jordan  (1960).  The  remainder  of  this  chapter  will  deal  with  the 
diffuse-hemispheric  irradiance  correlations  determined  from  the  1977 
through  1980  data  base  for  Fort  Collins,  Colorado. 

1.  Method  of  correlation 

The  daily  0IF/ETR  and  HEM/ETR  ratios  were  calculated.  For  each 
HEM/ETR  interval  of  0.1G  the  corresponding  values  of  DIF/ETR  were  ave¬ 
raged  and  these  mean  values  were  plotted  against  the  value  of  HEM/ETR 
for  the  mid-point  of  that  interval.  The  only  exception  to  this  approach 
was  for  the  HEM/ETR  interval  0.00-0.09.  Only  four  days  fell  into  this 
interval  with  values  of  0.07-0.09.  Therefore,  the  mean  DIF/ETR  value 
(0.08)  was  plotted  against  the  mean  HEM/ETR  value  (0.08).  A  least 
squares  polynomial  fit  was  then  obtained  for  these  means. 

The  cumulative  density  function  of  DIF/ETR  values  within  each 
HEM/ETR  interval  was  also  obtained,  and  the  20  and  80  percent  levels 
determined.  The  values  of  each  of  these  levels  were  again  plotted 
against  the  value  of  HEM/ETR  for  the  mid-point  of  that  interval  (except 
that  values  for  the  lowest  interval  were  again  plotted  against  a  HEM/ETR 
value  of  0.08),  and  a  least  squares  ploynomial  fit  obtained  for  each  of 
these  levels. 

2.  Correlation  equations 

The  interval  0.35  S  (HEM/ETR)  ^  0.80  included  approximately  86.3 
percent  of  the  total  number  of  days.  The  following  correlations  were 
obtained: 
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Mean 

-0.348  +  3. 93(HEM/ETR) 

-7.51(HEM/ETR)2  +  4. 14(HEM/ETR)3 

[3.2] 

20%  level 

-0.129  +  2.72(HEM/ETR) 

-5.8KHEM/ETR)2  +  3.42(HEM/ETR)3 

[3.3] 

80%  level 

0.243  +  0. 386(HEM/ETR) 

-0.124IHEM/ETR)2  -  0. 750(HEM/ETR)3 

[3.4] 

There  were 

no  points  in  the  interval  (HEM/ETR)  >  0.81. 

The  remaining  days  fell  in  the  interval  0.00  s  (HEM/ETR) 

<  0.35. 

correlations 

obtained  were: 

Mean 

0.010  -  0. 863(HEM/ETR) 

+  0. 200(HEM/ETR)2  -  1. 13(HEM/ETR)3 

[3.5] 

20%  level 

0.0913  -  0. 724(HEM/ETR) 

+  7. 62(HEM/ETR)2  -  12. 0(HEM/ETR)3 

[3.6] 

80%  level 

0.0006  +  1. 34(HEM/ETR) 

-  1.20(HEM/ETR)2  +  0. 273(HEM/ETR)3 

[3.7] 

It  interesting  to  note  that  while  32.4  percent  of  the  values  in  the 
Orgill  and  Hollands  (1976)  study  (based  on  four  years  of  data  for  To¬ 
ronto  Airport,  Canada)  fell  in  this  interval,  which  again  represents  ex¬ 
tremely  cloudy  conditions,  only  13.7  percent  of  Fort  Collins  days  were 
included.  Caution  should  be  used  in  any  direct  comparison,  however, 
since  their  work  was  based  on  hourly  values  and  the  present  study  deals 
only  with  daily  totals. 

Figure  2  shows  plots  of  the  above  correlation  equations.  A  compar¬ 
ison  of  the  correlation  equations  for  the  mean  values  in  the  present 
study  with  those  developed  by  Liu  and  Jordan  (1960),  based  on  ten  years 
of  data  for  Blue  Hill,  Massachusetts,  and  by  Orgill  and  Hollands  (1976), 


Figure  2.  Correlation  equations  between  daily  hemispheric/extraterrestrial  and 

diffuse/extraterrestrial  ratios  for.  means,  80  percent,  and  20  percent  levels. 
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based  on  four  years  of  data  for  Toronto  Airport,  Canada,  is  given  in 
Figure  3.  The  correlations  for  the  three  different  regions  are  very 
similar  for  HEM/ETR  values  less  than  0.15.  This  is  no  doubt  due  to  the 
fact  that  in  this  range  the  direct  component  is  zero  (or  very  nearly  so) 
and,  therefore,  the  total  insolation  is  very  nearly  all  diffuse.  In  the 
HEM/ETR  range  of  0.15  to  around  0.45,  the  present  study  gives  DIF/ETR 
values  that  lie  about  midway  between  the  other  two  correlations,  while 
above  HEM/ETR  values  of  0.60  the  present  study  and  that  of  Liu  and 
Jordan  (1960)  agree  quite  well. 


Figure  3.  Comparison  of  the  hemispheric/extraterrestrial  to  diffuse/extraterrestrial 
correlations  of  the  present  study  with  those  of  Liu  and  Jordan  (1960)  for 
Blue  Hill,  MA  and  Orgill  and  Hollands  (1976)  for  Toronto  Airport,  Canada. 


CHAPTER  IV 


SUMMARY  AND  CONCLUSIONS 

An  investigation  of  the  solar  radiation  environment  of  Fort 
Collins,  Colorado,  was  performed.  A  method  of  calculating  the  extrater¬ 
restrial  solar  radiation  over  various  time  scales  and  with  known  errors 

was  accomplished.  Uncertainty  of  the  extraterrestrial  values  resulting 

-2 

from  this  calculation  are  1.3  W  m  in  calculating  instantaneous  values, 

-2  -l  -2  -1 

4  W  hr  m  hr  for  hourly  values,  and  5  W  hr  m  day  for  daily  total 

calculations.  These  errors  are  certainly  less  than  those  to  be  expected 
in  the  measurement  of  surface  solar  radiation.  The  computer  program 
developed  to  perform  the  ETR  calculation  is  given  in  the  Appendix  and 
will  produce  daily  totals  for  one  year  for  about  $1.20  on  a  Control  Data 
Corporation  CYBER  system.  (Note:  This  is  not  intended  to  be  an  endorse¬ 
ment  of  the  CDC  CYBER  and  is  given  only  as  an  example  of  expected  cost.) 

Equations  correlating  the  diffuse  and  hemispheric  components  were 
developed.  Since  establishment  of  confidence  intervals  for  these  cor¬ 
relations  would  entail  assumptions  regarding  the  statistical  normalicy 
of  the  data,  correlations  for  the  20  percent  level  and  80  percent  level 
obtained  from  cumulative  density  functions  of  the  actual  data  were 
determined.  This  was  done  in  an  attempt  to  provide  the  maximum  infor¬ 
mation  with  the  fewest  assumptions.  Usinr  these  correlations  the  mean 
daily  diffuse  component  of  the  insolation  can  be  predicted,  given  the 
value  of  the  hemispheric  component.  These  present  correlations  were 
also  compared  to  those  found  by  Liu  and  Jordan  (1960)  and  Orgill  and 


Hollands  (1976)  and  were  found  to  agree  well  with  the  results  of  Liu  and 
Jordan  (1960)  except  in  the  HEM/ETR  range  of  0.20  to  0.60  where  the 
present  work  gives  values  of  DIF/ETR  ranging  from  0.01  to  0.04  higher. 
Except  in  the  very  lowest  values  ((HEM/ETR)  i  0.20)  the  correlations  of 
Orgill  and  Hollands  (1976)  consistently  over  estimate  the  DIF/ETR  values 
of  the  present  study  by  as  much  as  0.05.  While  this  study  was  based  on 
data  taken  only  at  Fort  Collins,  Colorado,  and  the  correlations  have  not 
been  applied  to  an  independent  data  set  for  other  locations,  the  author 
believes  the  correlations  given  should  prove  useful  for  other  climato¬ 
logical  ly-  simi lar,  high  altitude  sites. 

Finally,  an  investigation  into  extended  periods  of  low  insolation 
was  made.  For  the  four  year  data  set  six  periods  of  three  days  or  more 
were  found  where  the  diffuse  component  would  account  for  more  than  90% 
of  the  total  insolation  received.  Additionally,  during  a  58  day  period 
three  cases  totaling  29  days  were  found  in  which  the  direct  component  to 
extraterrestrial  direct  component  ratio  did  not  exceed  0.15.  The  four 
year  data  set  is  certainly  too  short  a  time  period  to  establish  a  clima¬ 
tology.  However,  the  results  of  this  study  tend  to  lend  creedence  to 
the  reputation  of  the  Colorado  Front  Range  as  being  very  suitable  for 
solar  energy  applications. 
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APPENDIX 


This  Appendix  is  a  listing  of  a  program  to  compute  daily  values 
of  extraterrestrial  solar  radiation  (HEMETR),  the  normal  component 
of  the  extraterrestrial  solar  radiation  (DIRETR),  and  the  length  of 
the  day  (DAYLEN). 
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PROGRAM  EXTRA 

COMPUTES  EXTRATERRESTRIAL  SOLAR  RADATKION  OK  BOTH  A  HORIZONTAL 
PLANE  (HF.METR)  AND  A  PLANE  NORMAL  TO  THE  SUN'S  RAYS  (DIRETR)  BY 
T.M.  THOMPSON'S  (1931)  METHOD  (MODIFIED  BY  J.R.  CON LEY ( 1982)) 
FOR  DETERMINING  SOLAR  POSITION.  (ETR  IS  IN  l /ATT  HOURS  PER 
SQUARE  METER).  ALSO  COMPUTES  DAYLENCTH  (DAYLEN). 

DAYLEN  UNITS  ARE  HUNDREDTHS  OF  HOURS  (IE.  932  IS  9.32  HOURS). 


INTEGER  ALLSOL 

DIMENSION  ALLSOL(3,31 ,12),  NDAY (12) 

COMMON  TIME I ,TIME2 ,SINLAT, COSLAT, XLONG, DECL,EQT,DIST, 

1  BEGGMT , EKDGI  IT, TIM INC , SOLC , ANGVEL, HEMETR , DIRETR , 

2  SUNR I S , S  UN  SET , D AYLEN , I YR , I MO , I DAY , I ET 
DATA  NDAY/ 31, 28, 31, 30, 31, 30, 3 1,3 1,30, 31, 30, 31/ 

0PEN( UNIT  =  7, FILE  =  'RADAT') 

INITIALIZE  ARRAY  ALLSOL  WHICH  WILL  CONTAIN  VALUES  OF 
EXTRATERRESTRIAL  SOLAR  RADIATRION  AND  DAYLENCTH. 


DO  1000  1-1,3 
DO  1000  J  =  1,31 
DO  1000  K  -  1,12 
1000  ALLSOL( I, J , K)  =  -9.00 

PROVIDE  LATITUDE  ,  LONGITUDE,  AND  STATION  NUMBER 
(EXAMPLE:  40  DEG  35  MIN  NORTH  =  40.5833  DEG,  40  DEG  SOUTH  = 
-40.0000,  90  DEG  WEST  =  90.0000  DEG,  90  DEG  EAST  = 

-90.0000  DEG). 


XLAT  =  40.5333 
XLONG  =  105.1333 
ISTA  =  53006 

CONVERT  LATITUDE  AND  LONGITUDE  TO  RADIAN  MEASURE 


PI  =  3.14159 
COKV  «  PI/ 180. 

XLAT  =  XLAT* CON V 
SINLAT  =  SIN(XLAT) 

COSLAT  =  COS (XLAT) 

XLONG  =  XLONC*CONV 

SET  SOLAR  CONSTANT  (WATTS  PER  SQUARE  METER)  AND  EARTH'S  ANGULAR 
VELOCITY  (RADIANS  PER  HOUR) 


C 


SOLC  =  1377. 
ANGVEL  =  15.*CONV 


o  o  o  n 
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SET  THE  TIME  PERIOD  (IN  DECIMAL  CUT  HOURS)  OVER  UHICtl  THE 
INTEGRATION  OF  ETR  IS  TO  BE  PERFORMED  (EX:  1010GMT  =  10.1667) 
AND  THE  YEAR  DESIRED. 

C 

BECCMT  =  11.00 
EtIDGMT  =  05.00 
IYR  =  7P> 

C 

C 

C 

C  SET  THE  INTERVAL  (IN  DECIMAL  HOURS)  OVER  M1IC11  INDIVIDUAL 
C  INTEGRATIONS  OF  ETR  ARE  TO  BE  PERFORMED. 

C  NOTE:  THIS  IS  THE  INTEGRATION  INTERVAL  (NOT  THE  PERIOD 
C  FOR  MUCH  A  FINAL  VALUE  IS  DESIRED)  AND  IS  NOR.'IALLY 

C  SET  TO  ONE  HOUR. 

C 

c 

TIM INC  »  1.0 
C 
C 
C 

DO  1500  IMO  =  1,12 
ND  =  NDAY(IMO) 

IF(MOD(  IYR, 4)  .EQ.  0  .AND.  IMO  .EQ.  2) HI)  =  2D 
DO  1500  I DAY  =  1 ,ND 
C 

C  THE  ADDITION  OF  0.5  TO  VALUES  BELOW  IS  DONE  TO  ROUND 
C  THEM  UP  BEFORE  FORCING  THEM  TO  BECOME  INTEGERS.  IN 

C  THE  CASK  OF  DAYLF.N  THE  MULTIPLICATION  FACTOR  OF  100 

C  FORCES  DAYLF.N  TO  BE  AN  INTEGER. 

C 

C 

ALLSOL( 1 , IDAY, IMO)  *  ETR ( )  +0.5 
ALLS0L(2,IDAY,  IMO)  =  HF.METR  +0.5 
ALLS0L(3, IDAY, IMO)  =  DAYLEN  *  100  +  0.5 
1500  CONTINUE 

WRITE(7 ,6000) ISTA, IYR 
DO  1600  J  =  1,31 

WRITE(7 ,6010)  J,((ALLSOL(I,J,K),K  =  1,12), I  =  1,3) 

1600  CONTINUE 
5000  CONTINUE 

6000  FORMAT( //// '  STATION  NO.',IC,5X,'  YEAR  19' ,I2//6X, 'DAY' ,6X, 

1  'JAN',7X, 'FEJ5',7X, 'MAR',7X, 'APR',7X, 'MAY' ,7X, 'JUN' ,7X, 

2  'JUL',7X,'AUG',7X,'SEP',7X, 'OCT' ,7X,'NOV' ,7X,'DEC') 

6010  FORMAT( 18 , 12110, 2( /<JX, 12110)) 

END 
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FUNCTION  ETR() 

THIS  FUNCTION  OPERATES  AS  THE  TIMINC  CONTROL  FOR  COMPUTATION  OF 
EXTERRESTRIAL  SOLAR  RAD I AT R ION  AND  TOTALS  DAILY  VALUES. 

COM! ION  TIME1 ,TIME2 .SINLAT, COSLAT, XLONG,DECL,EQT,DIST, 

1  BEGGMT.ENUCMT, TIMINC, SOLC.ANGVEL.HEMETR.DIRETR, 

2  SUNRIS,  SUNSET,  DAYLF.N ,  I YR ,  I  MO ,  I  DAY ,  I KT 


CALL  SUBROUTINE  TO  CALCULATE  ELAPSED  TIME  SINCE  BASE  DATE 
CALL  ELPSTM 


CALL  SUBROUTINE  TO  CALCULATE  SUN-EARTH  DISTANCE  PARAMETER 
CALL  DISTMN 


TIMEt  =  BECGMT 
HEMETR  =  0.0 
ETRDIR  =  0.0 
SUNRIS  =  0.0 
SUNSET  =  0.0 

IF(  ENDCMT  .LT.  BEGGMT)ENDGMT  =  ENDCtlT  +  24. 
10  CONTINUE 

IF(TIME1.GT.  ENDCMT)  GO  TO  20 
TIMF.2  =  TIME1+TIMINC 
HEMETR  =  ETR2()  +HEMETR 
ETRDIR  =  DIRETR  +  ETRDIR 
TIME1  =  TIME2 
GO  TO  10 
20  CONTINUE 

ETR  =  ETRDIR 

RETURN 

END 


n  rj  o  oonno  no 
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SUBROUTINE  ELPSTM 

COMPUTES  NUMBER  OF  DAYS  SINCE  THE  BASE  DATE  0  JAN  1974  FOR 
USE  IN  TIME  PARAMETERS  R1I0  AND  ALPHA. 

COMMON  TIMF.l  ,TIME2  ,  SINLAT  .COSLAT, XLONG,  DF.CL,  EQT,  DIST , 

1  BECGMT , ENDCNT, TIMINC , SOLC , ANCVEL , HEHETR ,DIRETR , 

2  SUNRIS, SUNSET, DAYLEN, IYR , I MO, I DAY, IET 


INITIALIZE  0  JAN  1974  TO  NUMBER  OF  DAYS  SINCE  1  JAN  1901 
I  BASE  =  26663 

COMPUTE  NUMBER  OF  DAYS  FROM  I  JAN  1901  TO  DATE  DESIRED 
LEAPYR  =  2 

IF( ( IYR/ 4*4)  .EQ.  IYR) LEAPYR  =  1 
IDYYR  =  IMO*275/9  +  IDAY  -  30 
IF ( IMO  ,G T.  2) IDYYR  =  IDYYR-LEAPYR 
IDYCT  =  ( IYR-1 )* 146 1/ 4  +  IDYYR 


COMPUTE  ELAPSED  TIME  IN  DAYS  FROM  BASE  DATE  TO  DESIRED  DATE 

IET  =  IDYCT-1BASE 

RETURN 

END 


o  o  n  o  on 
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SUBROUTINE  DISTMN 

COMPUTES  THE  SQUARE  OF  THE  RATIO  OF  THE  MEAN  TO  ACTUAL 
SUN- EARTH  DISTANCE 

COMMON  TIME  1 , TIME2 ,  SIN'LAT , COSLAT , XLONC ,  DECL,  EQT.DIST, 

1  BECCMT , ENDGMT , TIMINC , SOLC , ANGVEL , HEMET R, DI RETR , 

2  SUNRIS, SUNSET, DAYLEN, I YR, IMO, IDAY, IET 


COMPUTE  TIME  PARAMETER  ALPHA  FOR  LOCAL  NOON  OF  GIVEN  DAY 

XNOON  =  L2.+XLONC/ANCVE1 

ALPHA  *  -0172017221*(  IF.T+(XNOON/24.  )  ) 

5INAL  -  SIN (ALPHA) 

COSAL  -  COS( ALPHA) 

COS2AL  -  2 • *COSAL*COSAL- 1 « 

DIST  =  1.000145  +  0.033382*COSAL  +  0.001952*S INAL  + 

1  O.OOOA9F>*COS2AL 
RETURN 
END 


no  non  no  non  nooooo  no 


36 


FUNCTION  ETR2() 

COMPUTES  EXTERRESTIAL  SOLAR  RADIATION 


COMMON  TIME1.TIME2.SINLAT, COSLAT, XLONG, DECL, EQT, DIST , 

1  BEGCMT ,  ENDCMT ,  TIMINC ,  SOLC ,  ANGVEL ,  HEMETR , DIRETR , 

2  SUNRIS,SUNSET,DAYLEN, IYR,IMO, IDAY, IET 

OBTAIN  MIDPOINT  OF  CURRENT  TIME  INTERVAL  FOR  USE  IN  CALCULATING 
RHO 

TIME  =  (TIMEI  +  TI!!F,2)/2. 

COMPUTE  TIME  PARAMETER  RHO 

RHO  =  0.017 2027912*( IET  +  (TIME/24.)) 


COMPUTE  TRIGONOMETRIC  VARIABLES 
SINRHO  =  SIN (RHO) 

COSRHO  =  COS (RHO) 

SIN2RH  =  2.*SINRHO*COSRHO 
COS2RH  =  2.*COSRHO*COSRHO  -  1. 

SIN3RH  =  3.*SINRHO-4.*SINRHO*SINRHO*SINRHO 

COS3RH  =  4. *COSRHO*COSRHO*COSRHO  -  3.*COSRHO 

SIN4RH  =  4.*SINRHO*COSRHO  -  S.*SINRHO*SINRHO*SINRHO*COSRHO 

COS4RH  =  G.*COSRHO*COSRHO*COSRHO*COSRHO  -  S.*COSRHO*COSRHO  +1. 


COMPUTE  DECLINATION  AND  EQUATION  OF  TIME 

DECL  =  0.00666  -  0.4O0P.66*COSRHO  +  0.06523*SINRHO  -  0.00667* 

1  COS2RH  +  0.00057*SIN2RH  -  ().00274*COS3RH  +  0.00124* 

2  SIN3RH 

EQT  -  0.0090*COSRHO  -  0. 1227*SINRHO  -  0. 05 10*COS2RH  -  0.1581* 

1  SIN2RH  *  0.0056*SIN3RH  -  0.0021*COS4RH  -  0.0030*SIN4RH 


SINDEC  =  SIN (DECL) 

COSDEC  =  COS(DECL) 

HR 1  =  ANGVEL*(TIME1  -  12.  +  EQT)  -  XLONG 
HR2  -  ANGVEL* (TIHE2  -  12.  +  EQT)  -  XLONG 
HRO  =  ACOS(-(SINLAT*SINDEC)/(COSLAT*COSDF.C)) 

IF(HR1  .GT.  HRO)  GO  TO  10 

IF(HR2  .LT.  -HRO)  GO  TO  10 

IF(HR1  .LT.  -HRO)  HR1  =  -HRO 

IF(HR2  .GT.  HRO)  HR2  =  HRO 

ETR2  =  0.0 

DIRETR  =0.0 

IF(HR1  .GE.  HR2)  RETURN 

ETR2  =  S0LC*DIST*(SINLAT*SINDEC*(HR2-HR1)+ 

1  C0SLAT*C0SDEC*(SIN(HR2)-SIN(HR1 ) )) /ANGVEL 
DIRETR  =  SOLC*DIST*((HR2  -  HR1 )/ANGVEL) 

CALCULATE  SUNRISE  AND  SUNSET  HOUR  ANGLES  (DEFINITIONS:  SUNRISE  - 
CENTER  OF  SOLAR  DISK  IS  ON  THE  HORIZON  WITH  SUN  RISING;  SUNSET  - 
CENTER  OF  SOLAR  DISK  IS  ON  THE  HORIZON  WITH  SUN  GOING  DOWN). 

IF  (HR1  .EQ.  -HRO)  SUNRIS  =  HR1 
IF  (HR2  .EQ.  HRO)  SUNSET  =  HR2 

CALCULATE  DAY  LENGTH  FROM  SUNRISE  AND  SUNSET  HOUR  ANGLES,  ADDING 
CORRECTION  FOR  ATMOSPHERIC  REFRACTION  (-34  MINUTES  OF  ARC)  AND 
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THE  SUN'S  SEMIDIANETER  (-16  MINUTES  OF  ARC). 
DAYLEN  =  (SUNSET  -  SUNRIS)/ANCVEL  +  0.16 

RETURN 
10  CONTINUE 
ETR2  =  0.0 

DIRETR  =  0.0 
RETURN 
END 


